
C0(j3HHCKH yHHBepCHTeT „Cb. KjIHMCHT OxpHflCKH" 

Oh3hhcckh 4>aKyjiTeT 
KaTe^pa „ATOMHa 4)H3HKa" 



^HnjlOMHA PABOTA 
BjIIO'bK ^eTeKTOp OT CIJHHTHJiaipiOHHH 

BjiaKHa 3a HeyTpuHHaTa (J)a6pHKa 



PoceH HBaHOB MaTeB 



HayneH ptKOBO/piTeji: 
flon,. flc^H PyivieH I^eHOB 



Co4>hh, 2011 r. 



University of Sofia "St. Kliment Ohridski" 
Faculty of Physics 
Atomic Physics Department 



Master Thesis 

Scintillating fiber near detector for the 

Neutrino Factory 



Rosen Ivanov Matev 



Supervisor: 
Assoc. Prof. Roumen Tsenov, Dr. Sc. 



Sofia, 2011 



Abstract 



Neutrino Factory is a facility for future precision studies of neutrino 
oscillations. A so called near detector is essential for reaching the aimed 
precision of neutrino oscillation analysis. Main task of a near detector 
is to measure the flux of the neutrino beam. The brilliant neutrino 
source provides opportunity for precision studies of various neutrino 
interaction processes in a near detector. We present design concept 
of a scintillating fiber tracker capable of measuring the neutrino flux. 
The basic idea is to separate pure leptonic neutrino-electron scattering 
from the overwhelming background of neutrino- nucleon interactions. 
Monte Carlo simulation of the detector and simple event reconstruction 
algorithms are presented. It is shown that proper selection cuts increase 
signal to background ratio from ~ to 30-50 %. Two methods for 
extraction of signal events are developed. 
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Chapter 1 



Introduction 



Neutrino was postulated by W. Pauli in 1930 to explain the continuous kinetic energy- 
spectrum of electrons in beta decay. He claimed that the hypothetical particle was 
carrying away energy, momentum and angular momentum, so that conservation laws are 
held. The direct observation of the neutrino was not until 1956, when F. Reines and 
C. Cowan conducted their experiment at the Savannah River nuclear reactor in South 
Carolina. It was using neutrinos generated by the reactor and measured the 

z7 c + p-^e+ + n (1.1) 

interactions with their unique signature. 

1.1 Neutrino oscillations 
1.1.1 History 

The concept of neutrino oscillations was first stated in 1957 by B. Pontecorvo by analogy 
with the neutral kaons oscillations. The first observation of the phenomenon was made 
by R. Davis in the Homestake experiment in the late 1960s. The experiment counted the 
number of interactions of neutrinos emitted by nuclear fusion in the Solar core. Only 
one third of the theoretically predicted interaction rate was observed. Many subsequent 
experiments confirmed the deficit. In 2001, The Sudbury Neutrino Observatory provided 
the first direct evidence of solar neutrino oscillation. The experiment was able to detect 
solar neutrinos independently of flavor as well as only electron neutrinos. The measured 
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total flux of solar neutrinos agrees well with theoretical predictions, while electron 
neutrinos accounted for about 35 % of the total. 



1.1.2 Phenomenology 



To obtain a phenomenological description of neutrino oscillations one should abandon 
the Standard Model postulate of zero neutrino masses. In that case, neutrino mass 
eigenstates do not necessarily have equal eigenvalues and can in general not coincide 
with flavor eigenstates. Therefore, a flavor eigenstate of neutrino, \v a ) (I = e,/i, r), is in 
general expressed as superposition of mass eigenstates, (i = 1,2,3): 



(1.2) 



where U ai is an element of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) unitary 
matrix U. It is commonly parameterized with three mixing angles, 12 , 6> 23 and 6> 13 , and 
one complex phase 5: 



U 



li oW 
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C13 Si 3 e^ 
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^-siae^ ci 3 y 



(1.3) 



where c^- = cos% and = sin%. The evolution of an initial flavor eigenstate |i/ a ) 
according to quantum mechanics is given by the superposition of the evolved mass 
eigenstates: 



(1.4) 



where and pj are the energy and momentum of z/j, and L is the traveled distance. Due 
to the smallness of the neutrino mass, we can make the following approximations: 



t 



Ei = xlpj + mf « pi + 



m, 2 



(1.5) 
(1.6) 



where the natural units system (/i = c = 1) is used. The three mass eigenstates z/j have 
the same momentum because v a is produced with a definite momentum p m E. The 
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probability that vp is observed after a distance L is given by 

P(v a ^vp) = \(vp\v(t))? = 

= 5 a p - 4 MKiUajU^jUpi) sin 2 

i>j 

± 2 MKiUajU;^) sin I J , (1.7) 

i>j ^ ' 

where Am| = m 2 — m 2 is the mass squared difference between and Uj, and "+" is for 
neutrinos and "-" is for antineutrinos. The independent square mass differences are only 
two because of the obvious relation Amf 2 + Am 23 + Am 2 3 = 0. Therefore, the parameters 
describing neutrino oscillations in general are six: three mixing angles, one complex phase 
and two squared mass differences. Experiments have observed oscillations in two distinct 
regions, one at approximately 30 times larger E/L than the other. Therefore, one mass 
splitting is much smaller than the other two: Am 2 2 <C Am 2 3 ~ Am^. Because of this 
(and the smallness of sin# 13 ), oscillations in the regions E/L ~ Am 2 2 and E/L ~ Am 23 
can be approximated by oscillation of two neutrinos with effective mixing angle 9 and 
squared mass difference Am. The oscillation probability is then given by (for a ^ 0) 

P{v a ^vp) = sin 2 29 sin 2 • (1.8) 



E. m 4 L 

i 



AmfjL 
AE 



1.2 Neutrino Factory 

Neutrino Factory is a proposed facility for precise measurement of neutrino oscillation 
parameters. The Neutrino Factory will create an intense neutrino beam via decays of 
muons stored in racetrack-shaped storage rings. The decays of negative and positive 
muons 

jjT ->■ e~ + + V e (1.9) 
and /i + -> e + + + v e , (1.10) 

provide two separate neutrino beams. The baseline layout of the accelerator complex 
is shown on Figure [TTT] 1 1 1 . Two detector locations are currently envisioned, one at 
"intermediate" (3000-5000 km) baseline distance and the other at "long"( 7000-8000 km) 
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Neutrino Beam 




Muon Decay Ring ids-nf Boseime 2010/2.0 

Figure 1.1: Baseline layout of the Neutrino Factory accelerator complex. 

baseline distance. Specifications [I] give a total of 10 21 useful muon decays for both signs 
and both detector locations for 10 7 s3 useful time. 

1.2.1 Accelerator complex 

For achieving the baseline muon decays per year, a 4 MW primary proton beam is 
needed. The unprecedented beam power will be obtained by reusing and upgrading 
existing facilities and building new accelerators (linear or circular). The next step is 
to generate a secondary pion beam from interactions of the proton beam with a target. 
In order to avoid reinteractions of pions in the target, it has to be long and thin. This 
geometry imposes a great technical challenge due to the huge energy density created in 
the target. The solution being considered is a liquid jet mercury target |2]. Generated 
pions will be captured by a 20 T solenoid, which tapers down to 1.5 T in a distance of 
15 m. Contrary to a magnetic horn, the solenoid allows for both pion signs to be captured. 
The pion beam decays in a long decay channel, thus generating a tertiary muon beam 
with a large energy spread. The muons are bunched and phase-rotated in a sequence 

1 When "year" or "annual" is used further on, it refers to 10 7 s 
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of RF cavities so that a train of bunches with equal energies is formed. Following is a 
cooling channel, which reduces the transverse emittance of the muon bunches so that 
more muons fit in the acceptance of the following accelerators. Traditional techniques, 
such as stochastic, electron and laser cooling, are not applicable at the Neutrino Factory 
since they take in the order of seconds to minutes to cool a beam. Muons with a rest 
frame lifetime of 2.2 us need to be cooled with a method that is fast enough. A novel 
approach is considered called ionization cooling. The idea behind the method is to use 
ionization losses to reduce emittance in both longitudinal and transverse directions and 
then accelerate to restore the initial momentum. A dedicated experiment called MICE 1 3 
to prove the feasibility of the concept is in operation. Acceleration is divided into three 
stages, each chosen to use the most efficient acceleration method at its energy range. 
First, the muon beam is accelerated to 0.9 GeV energy by a linear accelerator. This is 
followed by two recirculating linear accelerators (RLA), which increase the energy of 
the beam to 3.6 GeV and 12.6 GeV, respectively. Finally the beam is accelerated in a 
fixed-field alternating-gradient (FFAG) accelerator. In both RLA and FFAG accelerators 
the beam makes many passes through the RF cavities, thus making them cost-effective 
solutions. The beam with an energy of 25 GeV is injected in two racetrack-shaped decay 
rings. They allow simultaneous storage of both muon signs. Each decay ring points 
toward one of the far detector locations. 

1.2.2 Physics 

The Neutrino Factory is the most promising future facility for studying of neutrino 
oscillation parameters. It has the following main goals: 

• accurately measure sin 2 (26*13) if it is relatively large, and have the potential to test 
very small values otherwise]^] 

• determine the hierarchy of neutrino mass eigenstates; 

• measure the value of the CP violating phase 5. 

The advantages of the Neutrino Factory compared to other present and future facilities are 
the intense flux, the ability to produce beams of both v (up) and V e (u e ), the cleanliness 
of the beams. Twelve oscillation channels are available for studies. The six channels 



2 Recently, the T2K experiment has shown indication that sin 2 (26*13) is large |4 . Best fits for normal 
and inverse hierarchies are 0.11 and 0.14, respectively. Statistical significance of the result is 2.5a. 
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available in the (u^, i/ e )-beam are shown in Table The other six are their charge 
conjugates. 

Table 1.1: Oscillation channels for a [y , z7 e )-beam decay mode). Oscillation channels in 
the other decay mode are charge conjugates of these. 



Channel 







appearance "Platinum channel" 






disappearance 


u , 




appearance 


v c 




disappearance 


v c 




appearance "Golden channel" 


v c 


—> V T 


appearance "Silver channel" 



The Golden channel u e — y u has a unique experimental signature of a muon with 
sign opposite to the muon signs in the decay ring. The channel gives opportunity for 
measuring both # 13 and 5 parameters. However, there is an intrinsic degeneracy that 
continuously many solutions to the equation P a p{9i3,S) = P exist. By studying both 
neutrino and anti-neutrino oscillations the continuous degeneracy can be eliminated. 
However, "intrinsic clone" degenerate solutions still remain. They arise from the lack of 
knowledge of the mass hierarchy and the octant where # 13 is located. These degeneracies 
can be resolved at the Neutrino Factory by studying different oscillation channels at two 
different baselines. 

1.2.3 Near detector 

A future neutrino facility will need near detectors in order to perform oscillation measure- 
ments with the required sensitivity. The near detector measurements that are essentialfor 
the neutrino oscillation analysis are 1 1 1 : 

• measurement of neutrino flux; 

• measurement of neutrino beam properties needed for the flux to be extrapolated to 
the far detector; 
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• measurement of charm production cross sections (charm production in far detectors 
is one of the principal backgrounds to the oscillation signal). 

The brilliant neutrino beam allows for other standard neutrino physics studies, such as 
measurement of cross sections, structure functions and sin 2 9w at neutrino energies in 
the 0-25 GeV range. 



Chapter 2 



Neutrino-electron scattering 



Neutrino-electron interactions involve only fundamental particles, thus they do not 
depend on unknown form-factors. Cross sections of these processes are straightforwardly 
calculated in the Standard Model framework. Any small uncertainties could come only 
from well measured Standard Model parameters. Therefore, such processes are suitable 
for measurement of neutrino beam fluxes, provided that beams are intense enough. The 
cross sections given in this chapter are calculated in the low-energy approximation with 
four-fermion interaction. This approximation is justified for interactions with s <C m\j. 
For a Neutrino Factory beam having maximum energy of E u = 25 GeV, the square of 
the invariant mass is 



Neutrino-electron interactions relevant for a Neutrino Factory beam are briefly described 
in the following sections. 

2.1 Quasi-elastic scattering 

The two purely leptonic interactions of neutrinos from the Neutrino Factory beam 
producing a muon in final state are 



s = 



(p* + Pc) 2 = 2E u m e + m 2 e & 0.03 GeV 2 . 




and z/ e + e~ — > ji~ + v p 



(2.2) 
(2.3) 
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The process in equation (2.2), known as inverse muon decay, has an isotropic differential 



cross section in the center of mass frame. The total cross section is given by (5] 

Gl (s -ml) 2 
a=^- (2.4) 

7T S 

where Gp is the Fermi coupling constant, s is the square of the invariant mass of the 



system and is the mass of the muon. The process in equation (2.3), known as muon 
production through annihilation, has differential cross section in the center of mass frame 
given by [5] 

= _<J {s-^f (£ , + ^ e) ^ + ^ 9) (2 5) 

and the total cross section is 

71 S 2 \ 3 



° = ~ ,2 E * E I> + O^Eu, ) ■ (2.6) 



The neutrino energy threshold (for electrons at rest) for both processes is 

2 2 
fjl^ — TYl 

E v > - « 10.9 GeV. (2.7) 

2m e 

The two quasi-elastic neutrino-electron scattering processes will be referred to with IMD. 



2.2 Elastic scattering 



There are two pure neutral current reactions of interest: 



^ + e ->^ + e" (2.8) 
and V/j, + e~ — > + e~. (2.9) 



Their respective cross sections are |6 



G F m e E u ( 9 16 A \ , 

a(u e-^u e) = — — 1 - 4 sin^ 9 W + — sin 4 9 W ] and (2.10) 

2n \ 3 J 

(- - ^ C|m c ^/1 4 . 2 16 . 4 _ \ 

er(z^e — >■ i^e) = - - - sin v w + — sin 6V I , (2.11) 



2tt \3 3 3 
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where 6\y is the Weinberg angle and E u ^> m c approximation was made. The other 
processes of interest result from both W and Z°-boson exchange: 

u e + e~^u e + e- (2.12) 

and V c + e~ — > V e + e~. (2-13) 

Their respective cross sections are |6 



a {v e e -» v c e) = G ^ Ev fi + 4 sin 2 6 W + ^ sin 4 9 W ) and (2.11) 
In \ 6 

v G\m c E v (\ 4 2/1 16 4/1 . 
cr(z/ c e — >■ i/ e e) = — - + - sin 6*^/ + — sin 6V . (2-15) 



2vr V 3 3 3 
The elastic neutrino-electron scattering processes will be referred to with ES 



2.3 Experimental signatures 



Total cross sections for neutrino-electron interactions as a function of neutrino energy are 



given on Figure 2.1 Despite their smallness, a couple of tons detector placed 100 m after 



the straight section end can provide sufficient interaction rate - Figure 2.2| However, 
inclusive charged current (CC) and neutral current (NC) neutrino interactions with nuclei 



u e + N^£ + X 
and v t + N ->• v t + X 



(2.16) 
(2.17) 



are a few orders of magnitude more frequent. For example, muon neutrino and anti- 
neutrino CC total cross sections in DIS regime are (7j cr/E u ^ fa 0.68 x 10 -38 cm 2 /GeV 
and a/E v ~ 0.33 x 10~ 38 cm 2 /GeV respectively. At 15 GeV the muon neutrino CC total 
cross section is ~ 1 x 10~ 37 cm 2 , compared to ~ 2 x 10~ 41 cm 2 for inverse muon decay 
— > fJL~v e . An obvious distinction between purely leptonic processes and processes 



(2.16) and (2.17) is the lack of hadronic system X in the former. Measurement of the 



recoil energy of the hadronic system can be used as good criterion for background 
suppression. Muons from quasi-elastic neutrino-electron scattering have a distribution 
peaked at very forward direction. At the Neutrino Factory, the polar angle of IMD 
muons does not exceed 5 mrad. The angular spread comes mainly from the intrinsic 
scattering angle ~4 mrad in processes (2.2) and (2.3), while neutrino beam divergence 
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Figure 2.1: Total cross section divided by neutrino energy for the leptonic interactions in a 
Neutrino Factory beam. Blue and green show processes which occur with fi~ in 
the storage ring, while red shows processes which occur with fi + in storage ring. 
Hexagon markers show CC interactions, triangle markers show NC interactions and 
rhombus markers show mixed CC+NC interactions. The threshold for quasi-elastic 
processes (blue) is evident at ~11 GeV. 



(and solid angle covered by detector) has little contribution. This kinematic property 
can be used as another event selection criterion. Polar angle distribution of electrons 
from neutrino-electron elastic scattering is ten times wider and is not suitable for event 
selection. On the other hand, the Q\Eg variablq^j provides good separation between signal 
and background for all neutrino-electron scattering processes. 



1 The variable 9fEg is proportional to Bjorken's y = 1 — Ei/E u (event inelasticity) for elastic and 
approximately proportional for quasi-elastic scattering. 
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Neutrino energy, E u [GeV] 



Figure 2.2: Number of neutrino-electron interactions for a nominal year of Neutrino Factory 
operation. Rates are calculated for 2.7 t detector with 1.5 x 1.5 m 2 frontal cross 
section and average Z/A ~ 0.54. Detector is placed 100 m after the straight 
section end. Dashed vertical line indicates threshold for quasi-elastic scattering. 



Chapter 3 
Near detector 



3.1 Measurement of neutrino-electron scattering 

To perform measurements of neutrino-electron scattering, the detector must be able 
to distinguish between the leptonic interactions and inclusive CC and NC neutrino 



interactions (2.16) and (2.17), which are a few orders of magnitude more frequent. To 



accomplish this task, the detector must make use of the distinctive event signatures 



mentioned in Section 2.3 This imposes certain general requirements to the design of the 



near detector. It should be/have: 

• massive enough to provide sufficient interaction rate; 

• tracker-like design to measure precisely the primary lepton's scattering angle; 

• magnetic field to measure the momentum of the primary lepton; 

• precise calorimeter to separate background events based on hadronic recoil energy; 

• low Z material to minimize multiple scattering and electromagnetic showering. 



3.2 Scintillating fiber tracker 



A detector design which fulfills the requirements stated in Section |3.1| is proposed. The 
detector design is essentially a scintillating fiber tracker with an incorporated calorimeter. 
A schematic drawing of the detectoiQ is shown on Figure 3^ A right-handed coordinate 



1 From now on, the scintillating fiber tracker will be referred to as "the detector". 
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20 modules 




-11 m 



5 cm ■ 



-7.2 mm 



Figure 3.1: Schematic drawing of the detector. 



system is adopted in which the origin is at the most upstream surface of the detector, the 
positive Z axis is aligned with the neutrino beam axis and the positive Y axis is pointing 
upwards. The detector consists of 20 square shaped modules placed perpendicular to 
the beam axis. Each module has a calorimeter section and a tracker section (also called 
tracker station). Modules are positioned equidistantly forming gaps filled with air. With 
larger distance between tracker stations, transverse displacement of hits is increased and 
thus angular resolution improved. The sides of the air gaps are covered with a layer of 
plastic scintillating bars. These layers are referred to as side slabs. The detector is placed 
in 0.5 T dipole magnetic field. Each station consists of one layer of fibers with horizontal 
orientation and another with vertical orientation. Each layer has four planes made of 
1 mm round fibers. They form a hexagonal pattern in the layer, thus minimizing dead 
volume. Calorimeter sections consist of plastic scintillating bars arranged in 5 slabs. All 
bars are oriented vertically, so that their axes are perpendicular to the magnetic field 
vector. Bars are co-extruded with a wavelength shifting (WLS) fibers inside and have 
10 mm by 30 mm cross-section. Both tracker fibers and WLS fibers in bars are read from 
both ends by silicon photomultipliers (SiPMs). Overall dimensions of the detector are 
~ 1.5 m x 1.5 m x 11 m. The detector mass is ~ 2.7 t. 



3.3 Silicon photon detectors 



SiPM photo detectors combine many small avalanche photodiodes operated in Geiger 
mode to form a single device (see Figure 3.2). They were proposed in the 1990s (8j|9). 
Detailed information about SiPMs can be found in a recent review paper |9]. In the past 
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Figure 3.2: Photograph of 1.2 mm diameter SiPM manufactured in FBK, Trento, Italy. 

several years the technology has undergone rapid development, making it competitive 
to conventional photomultiplier tubes (PMT). Advantages of SiPMs over PMTs are 
numerous and include: independence of external magnetic fields, high photon detection 
efficiency, low operating voltage, compact size. The most important advantage is that 
SiPMs can be manufactured in standard microelectronics facilities. Furthermore, it is 
possible to integrate electronics into the SiPM itself, lowering the cost per channel. This 
makes silicon photomultipliers especially suitable for application in detectors having 
hundreds of thousands fibers, such as the scintillating fiber tracker. A SiPM device is 
build from array of avalanche photodiodes (pixels). Size of pixels varies from 20 um to 
100 um. Each avalanche photodiode is operated in Geiger mode, meaning that a primary 
trigger (photon, thermal electron) creates a discharge in the pixel. This results to a gain 
in the order of 10 6 . Each pixel is rearmed by a resistor which quenches the discharges by 
the voltage drop induced from the discharge current. All pixels are connected in parallel, 
thus forming an analog signal which is a sum of all pixel signals. Therefore, the SiPM 
signal is proportional to the number of primary triggers, provided that it is not saturated. 
Although the SiPM acts as an analog device, in the case of few incident photons, the 
signal exhibits the discrete structure of the device, because all pixels are nearly the same. 
The operating voltage is in the 25-70 V range, a few volts above the diode's breakdown 
voltage. 



Near detector 



16 



The photon detection efficiency (PDE) of a SiPM depends on three factors: 

PDE = QE x e pt x e geom , (3.1) 

where QE is the wavelength dependent quantum efficiency, e pt is the probability of 
initiating a Geiger discharge by a primary trigger and e geom is the sensitive fraction of 
the device area. The gain of SiPMs is voltage dependent and varies across different 
production sets of the same device. Moreover, gain has also a temperature dependence, 
which additionally complicates calibration procedures. Discharges induced by thermal 
electrons are called dark counts. Dark count rate is temperature dependent and has a 
typical value of less than 1 MHz for 1 mm 2 devices at room temperature. When single 
photon counting is not essential, imposing a threshold on the signal can eliminate most 
of the dark counts. 



Chapter 4 

Monte Carlo simulation 



4.1 Neutrino beam simulation 

Neutrino flux at the near detector is generated by a Monte Carlo simulation of muon 



decays along the straight section of the NF decay ring 10 11 . This Section will briefly 
discuss the simulation. 

For the decay of the negative muon, the expected decay distribution for the Standard 
Model values of Michel parameters is |12| 



d?T 

~x 2 ((3-2x) +P M cos0(l-2x)) faru, (4.1) 



dxd cos 9 

Q -x 2 ((l-x) + P ll cos6(l~x)) foru e , (4.2) 
dxd cos U 

where x = 2E l/ /m II , P M is the muon polarization and 9 is the angle between the muon 
polarization vector and the direction of the neutrino. For the decay of the positive 
muon, the sign of the term containing cos# should be inverted. Energy and angular 
(tan -1 (dx / dz)) distributions of muons in the beam are taken to be Gaussian. It is 
assumed that the transverse distribution of muons in the beam pipe is a Dirac delta 



function. The parameters used for the beam simulation are given in Table 4.1 



The simulation of each muon decay goes as follows. First, a position of the decay 
along the straight section is drawn from a uniform distribution. Then, for each neutrino, 



x and cos# is drawn from distributions 4.1, 4.2 using the acceptance-rejection method. 



An azimuthal angle <fi is drawn from uniform distribution. Having x, cos 6 and (f>, the 
momentum of the neutrino in the center of mass (CM) frame of the muon is obtained. 
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Table 4.1: Parameters of the beam in the muon storage ring used in the simulation. 



Parameter 


Value 


Muon energy central value 


25 GeV 


Muon energy RMS 


80 MeV 


Muon angular spread 


0.5 mrad 


Muon polarization 





Straight section length 


600 m 




'0 5 10 15 20 25 ' 10 20 30 40 50 60 70 

E v [GeV] Distance from axis [cm] 

Figure 4.1: Energy spectrum (left) and radial distribution (right) of neutrinos on plane per- 
pendicular to the straight section and 100 m away from its end. The distributions 
are normalized to the baseline specification of 2.5 x 10 20 muon decays per muon 
charge per straight section. Dashed line indicates the threshold for the leptonic 
processes with muon in the final state. 

Energy and polar angle of the muon are drawn from the according Gaussian distributions. 
Combining with an uniform random azimuthal angle, one obtains the muon momentum 
vector in laboratory frame. The neutrino momenta vectors in the CM frame are boosted 
to obtain the momenta in the laboratory frame. Electrons (positrons) from the decay 
are not simulated. 



The most upstream plane of the near detector is to be situated at 100 m after the 
straight section end. The energy and space distributions of neutrinos at that plane is 
shown on Figure 4.1 
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4.2 Neutrino interactions 



The neutrino interactions with nuclear targets and electrons is simulated with the GENIE 
neutrino Monte Carlo generator |l3j. The following neutrino interaction processes are 
simulated by GENIE: 

• Quasi-elastic scattering; 

• Elastic NC scattering; 

• Baryon resonance production in charged and neutral current interactions; 

• Coherent neutrino-nucleus scattering; 

• Non-resonant deep inelastic scattering; 

• Quasi-elastic charm production; 

• Deep-inelastic charm production; 

• Neutrino-electron elastic scattering; and 

• Inverse muon decay. 

The process z7 e + e~ — > + jjT is not included in GENIE. For the simulated beam, the 
process above gives 4 times less interactions than inverse muon decay as can be seen 



from Figure [272} For this work GENIE version 2.6.2 has been used. Configuration was 
changed, so that short-lived charm particles are decayed by GENIE. Apart from this 
change, default values were used for all other parameters. 

The NF neutrino flux is fed to GENIE by implementing a flux driver interface. The 
detector description used for event generation is a simplified version of the actual detector 
geometry. Each detector module is implemented as a solid box made of polystyrene. 
Air gaps are left as they are. Covering side bars are not implemented as interactions 
there are not in fiducial volume. For each neutrino interaction, GENIE produces a list 
of intermediate and final state particles with their kinematic characteristics. The type 
of the process involved in the interaction as well as some kinematic variables are also 
provided. 
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4.3 Detector simulation 



4.3.1 Particle transport 



For the simulation of detector response to neutrino interactions, the Geant4 14,15 
software platform was used. Geant4 is a toolkit for "the simulation of the passage of 
particles through matter". From the user's point of view, to simulate an event with 
Geant4 one has to define (by implementing the toolkit's interface classes): 

• primary particle generator; 

• list of particles to be transported (including decay products and particles produced 
in interactions) through the detector and relevant physics processes (physics list); 

• detector description in terms of geometric volumes and comprising materials; 

• "sensitive" volumes, where "snapshots" of the physical interaction and state of 
particles are recorded to produce "hits". 

All final state particles produced by GENIE are passed to Geant4 simulation as 
primary particles. The physics list used in the simulations is "Simple and Fast Physics 
List" fl6| . The scintillating fiber tracker geometry was implemented in a flexible manner, 
so that detector parameters are set with a configuration file. Plastic scintillating fibers 
and bars are defined as sensitive volumes. The hits produced in sensitive volumes contain: 
unique id of the volume, id of the particle, energy deposition, time and position of the 
interaction. True particle trajectories are recorded for some of the particles. Part of the 
information recorded in the simulation is used to produce electronic signals (digitization). 
Other data is retained for later comparison of reconstructed data with true data (often 
called "Monte Carlo" truth). By doing the comparison with MC truth, the performance 
of the detector design and reconstruction algorithms can be analyzed. A basic event 
display has been developed to enable visual inspection of events. Part of an example 



event display is shown on Figure 42 Apart from MC truth, the event display shows 
some reconstructed data. 

Tracking of individual particles with Geant4 is stopped when their energy is below 
a certain threshold and it is guaranteed that they will not cross a volume border. 
Scintillation process was not included in the Geant4 simulation because tracking of 
individual scintillating (optical) photons is a very time consuming task. 
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Figure 4.2: Example event display showing a CC resonance production event. 

4.3.2 Digitization 

Digitization algorithms take the output hits of the Geant4 simulation and transform 
them into electronic signals. 

Fiber signal digitization 



As simulation of optical photons is not feasible in event by event basis, a toy Monte Carlo 
algorithm was developed to simulate the processes of scintillation, photon propagation 
into fibers and SiPM response. For each energy deposit in a fiber, a random number of 
scintillating photons is drawn according to the scintillation yield of the material. The 
scintillating photons are emitted isotropically, hence not all of them suffer total internal 
reflection from the walls of the fibers. The fraction of photons which do not exit the 
fiber is called trapping efficiency. It is determined from refraction indexes and fiber 
cross-section. For round fibers manufacturers provide a "minimal" trapping efficiency, 
namely the efficiency for photons propagating in a single central plane. However, it has 
been shown (T7j that the real trapping efficiency is about two times higher than that. A 
value of 8 % for the trapping efficiency is used in the simulation. The trapping efficiency 
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and light attenuation are taken into account to produce a number of photons arriving 
at the end of the fiber. If there are multiple hits in a fiber, photons arriving at the 
fiber end are summed up. The obtained number of photons is multiplied with an optical 
coupling efficiency factor to get the number of incident photons on the SiPM surface. 
Multiplying the last number by the photon detection efficiency (PDE) one obtains the 
number of primary avalanche triggers. Dark counts are simulated by further adding a 
Poisson distributed number to the primary triggers. The final electronic signal response 
is parameterized by a function that takes into account cross-talk effect and the single 
pixel response distribution. Only signals with amplitude larger than threshold equivalent 
to 2.5 fired pixels are recorded. For each fiber, the signals from both photo detectors 
and an id of the fiber are stored in a structure called fiber digit only if both signals are 
above the threshold. It seems dark counts are not a major problem for 1 mm fibers and 
modern SiPMs, which have less than 1 MHz/mm 2 dark count rate at room temperature. 
With the imposed signal threshold and coincidence requirement the number of "false" 
fiber digits is under 1 % of all fibers in the detector. Therefore, the probability to have 
two or more adjacent "false" fiber digits which mimic true particle hits is negligible. 

Bar signal digitization 

As above, Geant4 simulation does not treat optical photon processes in the scintillating 
bars and the WLS fibers inside. Ideally, one can make a detailed simulation or an experi- 
ment with prototype bar and parameterize the signal response (wavelength spectrum) 
at the end of the WLS fiber. However, we do not have such parameterization. For this 
work, bar signals were naively digitized by smearing energy deposits with 20 % Gaussian. 
Signals below a threshold of 0.5 MeV equivalent are discarded. For each bar, the signals 
from both photo detectors and an id of the bar are stored in a structure called bar digit. 



Chapter 5 



Software reconstruction 

Raw output from both Monte Carlo simulation and a real experiment is a set of digi- 
tized signal amplitudes. To understand what processes took place and what particles 
propagated through the detector, one needs to obtain measurements of relevant physical 
quantities characterizing them. The aim of reconstruction algorithms is to infer high 
level information about the event from the raw detector output. 

5.1 Initial processing 

The first step of reconstruction process is to obtain some basic physical quantities from 
the digital signals. First, signal from the photo detectors should be translated back to 
number of incident photons. For fibers in tracker stations this quantity is sufficient. On 
the other hand, for bars, one should further obtain an absolute measurement of the 
energy deposit via a dedicated calibration procedure. Since this is a simulation rather 
than real experiment, the above procedure is assumed perfect and the calibration step is 
skipped. 

To simplify further the reconstruction, neighboring fired fibers are grouped into 
clusters. Two fibers are considered neighbors if their projections on the XY plane are 
either neighboring or overlapping. Depending on the fibers orientation, (x, z) or (y, z) 
position of the cluster is calculated by taking weighted average of the fiber positions. 
Weighted standard deviations from the mean coordinates are also calculated. A cluster 
object then contains list of fibers, two coordinates with their standard deviations and an 
amplitude, which is the sum of the numbers of incident photons on the photo detectors 
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Figure 5.1: Vertex volume reconstruction example. The two adjacent most upstream stations 
with at least one cluster each are considered. The clusters (pink ellipses) with 
highest X coordinate are selected from both stations and a line passing through 
their centers is taken. Another line is constructed from the clusters with lowest X 
coordinates. The most upstream fired bar, which is between the lines (but not 
upstream of the intersection point) is identified as vertex bar (purple box). 



at the fibers. Clusters composed of only one fiber are discarded to eliminate false hits 
due to SiPM noise. 



5.2 Vertex reconstruction 

The point where a neutrino interaction takes place is called vertex. The scintillating 
fiber tracker, like a typical neutrino detector, serves as a target. Therefore, one does not 
have a priori knowledge of the precise vertex coordinates. The vertex reconstruction 
algorithm attempts to find the volume (bar or fiber), where the vertex is located. Only 
the two most upstream stations with at least one cluster are considered. Since particle 
tracks are not curved in the XZ plane, we can linearly back-project the positions of the 
clusters towards the event vertex. To do this, the clusters with highest X coordinate are 
selected from both stations and a line passing through their centers is taken. Another 
line is constructed from the clusters with lowest X coordinates. The most upstream fired 
bar, which is between the lines (but not upstream of the intersection point) is identified 



as vertex bar (Figure 5.1). If there is no such bar, the cluster with the highest amplitude 
in the first upstream station is picked and its most upstream fiber is labeled as vertex 
volume. If there are no two adjacent stations with at least one cluster in each, the vertex 
cannot be reconstructed and the event is discarded. 86 % of the events have their vertex 
volume reconstructed. 
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Figure 5.2: Vertex Z coordinate resolution. Distributions of the difference of reconstructed 
and true vertex Z coordinate are shown for CC u interactions (solid blue), CC v e 
interactions (dotted red) and NC {y„ , V e ) interactions (dashed green, x 5). Vertical 
dashed lines indicate Z boundaries of a bar centered at 0. Uniform distribution of 
events with correctly identified vertex bar is clearly visible. Events in the narrow 
peak centered at have their vertex volume correctly identified in a cluster. The 
shoulder on the right (and the dip left of 0) for CC V e interactions (red) correspond 
to events with a reconstructed vertex bar adjacent and downstream of the true 
vertex bar. This misidentification may be explained with energy deposition in 
the vertex bar below the threshold and is only visible in antineutrino events. 
Misidentification probabilities for the three samples are shown in the legend. 




However, the vertex volume is not always correctly identified. The performance of 



the algorithm is shown on Figure 5J2 in terms of difference between Z coordinate of the 
reconstructed volume's center and the true vertex Z coordinate. Difference distributions 
are shown for interactions of (v , I?e)-beani. 



5.3 Primary track reconstruction 

The scintillating fiber tracker is a very fine grained detector. Because of this, recon- 
struction of the full event topology might be attempted. However, for the purpose of 
measuring purely leptonic processes, only events with a single track originating from the 
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vertex are considered. The main goal of reconstruction is to identify the primary (lepton) 
track and measure its parameters. 

Only events with one cluster per orientation in the first two stations are considered. 
This selection has low efficiency, because an event with a single primary particle might 
still have more than one cluster in the either of the first stations. (Consider 5-electrons 
or interactions of bremsstrahlung photons) Furthermore, as both angle and momentum 
of the primary track need to be measured, only events with hits in at least three tracker 
stations are retained. The remaining fraction of leptonic events after the above selections 
is 61 % for IMD and 52 % for ES. Angle and curvature of tracks is measured using only 
hits in the first three stations. In third station, if there are more than one clusters, only 
the cluster nearest to the straight line defined by clusters in the first two stations is 
considered. At the end, for the primary track we have: initial position, initial slopes 
{x'q = dx/dz and y' = dy/dz), initial momentum and charge. 

To asses the performance of the primary track reconstruction algorithm, the recon- 
structed track initial angles and momenta are compared to MC truth. Distribution of 
difference between reconstructed initial angle of the primary track 9 rec and true initial 



angle of the primary lepton 9 true is shown on Figure 5.3 From the Gaussian fit to the 
distribution, it is seen that the reconstructed angle is unbiased and the resolution (cr 
parameter of the fit) is ~ 0.5 mrad for both muons and electrons. Distribution of relative 
difference between reconstructed initial momentum of the primary track and true initial 
momentum of the primary lepton 

APrec Ptrue / c- i \ 

= , (5.1) 

Ptrue 



is shown on Figure |5.4| For muons, the Gaussian fits show that the reconstructed 
momentum have only small bias (due to energy loss) and resolution is going up to 
~ 9 % for the highest energy muons. For electrons, the distribution is biased towards 
the negative values with a heavy negative tail. The reason for this underestimation 
of electrons' initial momenta is that they loose momentum due to bremsstrahlung and 
ionization. If energy loss is taken into account with Kalman filter fitting, for instance, 
bias can be reduced. Kalman filter fitting [18] is a powerful method, which can improve 
angular and momentum resolutions. This is true especially if full advantage of pattern 
recognition abilities of the method are used. 
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Figure 5.3: Obtained angular resolution for muons (left) and electrons (right). Gaussian fits 
are shown with dashed lines. The sample mean, standard deviation and the fit 
parameters are shown in upper left corners. 
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Figure 5.4: Obtained momentum resolution for muons (left) and electrons (right). Gaussian 
fits are shown with dashed lines. The sample mean, standard deviation and the 
fit parameters are shown in upper corners. For muons, distribution is shown for 
two samples of events: one with true muon momentum in the [10-12] GeV range 
(blue) and one in the [20-25] GeV range (green). 
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5.4 Spatial distribution of energy deposits 

Having this rather simple event reconstruction, a backward going particle will cause 
wrong reconstruction of vertex position. Moreover, the backward going track will be 
wrongly identified as a primary lepton track. Some background events have particles 
which do not cross the first tracker station after the vertex and do not deposit enough 
energy in the vertex volume. In order to select events with a single primary particle one 
might exploit the spatial symmetry of energy deposits relative to the track. A variable 
which has good signal/background event separation power is the following. Firstly, we 
calculate the energy deposit weighted average position of the hits in the bars in the i-th 
slab: 

obs _ Ej XjjAEjj 

where j denotes the bar in the slab, is the X coordinate of the center of the bar and 
AE^ is the energy deposit in the bar. Secondly, the expected position of the primary 
track at the i-th slab is calculated: 

xT P = %o + x' (zi - Zq), (5.3) 

where Zi is the Z coordinate of the i-th slab's center, xq and z are the primary track's 
initial coordinates and x' = dx/dz is the primary track's initial slope. Finally, the 
variable is 



%dev 

n 

i 

where the summation is over all fired slabs. Such variable is expected to have low values 
for events with a single muon, medium values for electron showers and high ones for high 
multiplicity asymmetric events. 

Overall reconstruction efficiency (remaining signal events) is 53 % for IMD and 45 % 
for ES. Inefficiencies come from inability to reconstruct the vertex, selection of events 
based on the number of clusters in the two most upstream fired stations and inability to 
reconstruct primary track's initial angle and momentum. 



Chapter 6 
Data analysis 



6.1 Neutrino-electron scattering event selection 



The Neutrino Factory beam is composed of v and u e neutrinos in fi decay mode and 
from z7 M and u e neutrinos in fi + decay mode. Since we can distinguish interactions only by 



their final state, we should group neutrino-electron interactions as follows: processes (2.2) 



and (2.3) in the IMD group; processes (2.8) and (2.13) in the EST group; processes 



(2.9) and (2.12) in the ES^ group. Sometimes ES will be used when referring to 



both ES - and ES + . Three neutrino-electron event samples are selected from all events 
corresponding to the three groups. IMD and ES~ samples are selected in fi~ decay mode, 
while ES + sample is selected in fi + decay mode. 



6.1.1 Calorimetric vertex selections 

Both IMD and ES events have a property of low (consistent with single particle) energy 
deposition near the vertex. To exploit that, a cut on energy deposit in the vertex bar of 
4 MeV is imposed. If vertex is in a cluster of fibers, a cut is made on its amplitude (the sum 
of fiber amplitudes). In some background events, energetic (hundreds of MeV) charged 
hadrons escape through the air gaps leaving small or no deposits in the calorimetric 
sections. Therefore, it is required that there is no activity in side bars covering air gaps 
adjacent to the vertex. Another vertex related cut is the requirement that there is no 
energy deposits upstream of the vertex (backward deposits). Distributions of energy 
deposit in the vertex bar and energy deposit in side slab for signal and background events 



are shown on Figure 6.1 and Figure 6.2 
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Figure 6.1: Distributions of energy deposit in vertex bar for IMD (blue), ES~ (red) and 
background (black) events in mode. The fraction of events contained in the 
plot is indicated in the legend. All distributions are normalized to a unit area. 
This plot shows only events in which the reconstructed vertex volume is a bar. 
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Energy deposit in side bars near the vertex [MeV] 

Figure 6.2: Distributions of total energy deposit in bars covering air gaps adjacent to the 
vertex for IMD (blue), ES~ (red) and background (black) events in [i~ mode. The 
fraction of events contained in the plot is indicated in the legend. All distributions 
are normalized to a unit area. 
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Figure 6.3: Distributions of (AE*) (left) and max(AEi) (right) for IMD (blue), ES~ (red) 
and background (black) events in fi~ mode. The fraction of events contained in 
the plot is indicated in the legend. All distributions are normalized to a unit area. 



6.1.2 Other calorimetric selections 

To select IMD events, one can rely on the specific properties of muon dE/dx. At energies 
in the range of 11-25 GeV, a muon is nearly a minimum ionizing particle. Therefore, the 
following cuts were applied when selecting IMD events: mean of all slab deposits (AJ5j) 
is less than 3 MeV and the maximum max(AEi) is less than 12 MeV. For ES events, 
primary electron might induce shower in the detector, thus such cuts are not suitable. 



Signal and background distributions of (AEi) and max(AEi) are shown on Figure 6.3 
For selecting IMD events, the x& ev variable is cut at 2.5 cm, while for ES events a more 
relaxed cut at 15 cm is applied. Signal and background distribution of the Xdev variable 
are shown on Figure |6.4| 



6.1.3 Kinematic selections 

As signal events always have negatively charged primary lepton (fi~ or e~), a cut q/p < 
on the primary track is imposed. Quasi-elastic neutrino electron scattering has a threshold 
at ~11 GeV. Therefore, when selecting the IMD sample, events with primary track 
momenta of less than 10 GeV are discarded. In order to avoid contamination of IMD 
events in the ES~ sample, only events with primary track momentum less than 10 GeV 
are retained in this sample]^] This is a reasonable cut, since majority of ES - events 

: A full containment calorimeter enclosing the detector will be able to distinguish electrons from muons 
with high confidence. However, such additional calorimeter is not envisaged in the current near 
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Figure 6.4: Distributions of the Xdev variable for IMD (blue), ES _ (red) and background 
(black) events in fj,~ mode. The fraction of events contained in the plot is indicated 
in the legend. All distributions are normalized to a unit area. 



have primary track momentum of less than 10 GeV. Distributions of primary track's 



reconstructed q/p and momentum are shown on Figure 6.5 



The overall efficiencies of the selection cuts are 89 % for IMD sample, 73 % for 
ES~ sample and 87 % for ES + sample. All cuts with their efficiency and puritj^] are 
summarized in Appendix [X} 



6.2 Background subtraction 



Absolutely clean sample of signal events cannot be selected with a reasonable efficiency 
by employing selection cuts only. Therefore extrapolations of certain background distribu- 
tions should be made in order to subtract background from event samples. Background 
subtraction is different from event selection in the sense that after it is done, one obtains 
a distribution of signal events and loses event by event information. 



detector design and the author prefers not to make assumptions on electron/muon discriminating 
capabilities. 
2 Purity is the ratio ^ W '!r°' . 
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Figure 6.5: Distributions of primary track reconstructed q/p (left) and momentum (right) for 
IMD (blue), ES - (red) and background (black) events in /i~ mode. The fraction 
of events contained in the plot is indicated in the legend. All distributions are 
normalized to a unit area. 



We have chosen to do background subtraction in terms of primary lepton (/i or e) 
kinematic variables for the following reasons: 

• they provide the most powerful signal/background separation criteria; 

• background distributions can be most reliably extrapolated in terms of kinematic 
variables; 

• we can reliably measure primary lepton angular and momentum resolutions. 

The primary lepton 's scattering angle Og. an d initial momentum p are measured in 
the detector. Primary muons with energy above 1 GeV and all primary electrons are 
ultrarelativistic. Therefore, it is justified to approximate the primary lepton's initial 
energy with E = pc. For simplicity, we will perform one dimensional analysis, where 
histograms over a variable x = f(9,p) will be used for background subtraction. In the 



case of IMD signal extraction, scattering angle 9g and 9fEg variable (see Section 2.3) can 
be used to discriminate signal from background. In the case of ES signal, background is 
well separated only when one expoloits 9\Eg variable. The distributions of 9g and OjEg for 



IMD signal, ES~ signal and background are shown on Figure |6.6| Background distribution 
over 9\Ei variable is nearly flat. This fact allows for a simple parameterization of the 
background distribution. Our further analysis is made using only 9\Eg variable. 

As a first approximation we assume that the energy distribution of neutrinos in the 
beam is known. Then, by measuring only the number of signal events, one can make 
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Figure 6.6: Distributions of reconstructed primary track initial angle 0£ and d\Eg variable for 
IMD (blue), ES~ (red) and background (black) events in fjT mode. The fraction 
of events contained in the plot is indicated in the legend. All distributions are 
normalized to a unit area. 



absolute flux normalization. Two methods of obtaining the number of signal events are 
discussed below. 



6.2.1 Linear fit method 

Linear fit method relies on the nearly flat shape of the respective background distribution. 
The idea is to estimate the background under the signal peak by linear extrapolation 
from signal-free region. First, an interval on the Q\Ei variable is defined by the following 
rules: 

• the interval low limit is close to the signal peak; 

• there are almost no signal events in the interval (according to MC simulation); 

• the background is approximately linear in the interval. 

Then, the histogram is fitted with a straight line in the interval. Finally, the line is 
extrapolated towards zero to estimate the number of background events under the signal 
peak. The histograms over Q\En and the linear fits for the three event samples under 



consideration are shown in Figures 6.7, 6.8 and 6.9 Comparison between the estimated 



and the true number of signal events is given on Table |6.1[ It is seen that the true 
values lie within the 95 % confidence intervals of the predictions. This means that the 
hypothesis that the predictions are unbiased cannot be excluded. The main source of 
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Figure 6.7: Distributions over O^E^ for the IMD sample. The leptonic events histogram is 
filled with solid gray, the hadronic events histogram is hatched and the total 
spectrum is in black. The two cuts bounding the fit interval are drawn with 
dashed line. The red line indicates the background extrapolation. 



systematic error related to the method is the assumption about the linearity of the 
background shape. From the fit of the ES~ and ES + samples it is estimated that the 
systematic error of the method 

I Af si 9 _ K[ si 9 I 

A I fit 'true I /. - \ 

syst = NStg , (6-1) 



is less than 1 %. However, to give conclusive estimation of the systematic error, one 
should investigate if and how various parameters of simulation and selections influence 
the background shape. 



6.2.2 //+-method 

IMD interactions are present only in the \l~ decay mode. The idea of the /i + -method 
is to estimate the background under the IMD signal peak exploiting the distribution of 
positive muons detected in (z7 M , z/ e )-beam. The idea was previously used for measuring 



IMD process in the CHARM II detector 19 20 . In the near detector, an event sample 
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Figure 6.8: Distributions over 9^E e for the ES _ sample. The leptonic events histogram is filled 
with solid gray, the hadronic events histogram is hatched and the total spectrum 
is in black. The two cuts bounding the fit interval are drawn with dashed line. 
The red line indicates the background extrapolation. 
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Figure 6.9: Distributions over 9^E e for the ES + sample. The leptonic events histogram is filled 
with solid gray, the hadronic events histogram is hatched and the total spectrum 
is in black. The two cuts bounding the fit interval are drawn with dashed line. 
The red line indicates the background extrapolation. 
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Figure 6.10: Left plot shows distributions over O^E^ for the IMD sample. The leptonic events 
histogram is filled with solid gray, the hadronic events histogram is hatched 
and the total spectrum is in black. The two cuts bounding the normalization 
interval A are drawn with dashed line. The red line indicates the normalized /i + 
histogram. Right plot shows ratio of the /j>~ histogram and the fi + histogram 
over O^Efj,. The errors are statistical and bars correspond to la. Horizontal 
dotted line indicates the normalization ratio. 



from the (z/ M , i/J-beam events is selected with the same selection cuts as the IMD 
sample. For consistency, we again consider the O^E^ distributions. Muon antineutrino 
inclusive CC interactions V^N have approximately twice higher cross section than v N 
interactions. To account for that, the O^E^ histogram for fi + should be normalized to 
the O^Ep histogram for jjT. First, the ratio of the \x~ background histogram and the fi + 



histogram is calculated - Figure 6.10 (right). It is seen, that the ratio under the signal 
peak (O^E^ < 0.4 rad 2 MeV) is at the same level as the ratio outside the signal peak. An 
interval outside the IMD signal peak and with approximately constant ratio of and 
yU + -events is defined: A = [0.4, 2.0] rad 2 MeV. The normalization factor is then 

N A 

R= -4- = 0.210 ±006, (6.2) 
where N A (N^ + ) is the number of events in the interval A for the (/i + ) histogram. 



The \i histogram and the normalized fi + histogram are shown in Figure 6.10 (left). Apart 



from normalization, there are second order differences between fi and fi + distributions 



19,20 , which arc not taken into account. 
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Table 6.1: Estimated number of signal events for the three event samples. The result in the 
last row was obtained using the fi + background subtraction method, while the 
other three results were obtained using linear fit background subtraction method. 
Statistics correspond to 2.3 x 10 19 n~ decays and 2.3 x 10 19 /i + decays, which is 
approximately a tenth of the nominal year. 



Event 


Selection 


Overall 


Purity 


All events 


Signal events 


Signal events 


sample 


eff. 


eff. 








from fit 


IMD 


86 % 


46% 


81 % 


3520 


2850 


2926 ± 59 


ES~ 


70 % 


32% 


61 % 


7355 


4491 


4479 ± 86 


ES+ 


83 % 


37% 


63 % 


16964 


10607 


10512 ± 131 


IMD 


86 % 


46% 


81 % 


3520 


2850 


2831 ± 61 



Chapter 7 



Conclusions 



We have proposed a conceptual design for the tracker part of the Neutrino Factory near 
detector. A full Monte Carlo simulation of neutrino interactions and particle transport 
was implemented using standard tools. An idealized detector geometry and simplified 
signal digitization was used in the MC simulation. Basic reconstruction algorithms were 
developed to extract information relevant for neutrino-electron scattering measurements. 
Selection procedures, which increase signal to background ratio from ~ 1CT 4 initially to 
~ 30-50 %, are defined. One dimensional signal extraction analysis was performed. It is 
demonstrated that the number of neutrino-electron scattering events can be measured 
from the 9jEe variable histogram. A direct comparison between measured and true 
number of signal events shows a deviation of no more than 1 %. It is worth noting 
that MC truth was not used in reconstruction and signal extraction. Therefore, the 
results obtained in this work are a reliable approximation to results obtainable at a real 
experiment, provided that MC simulation is close to reality. 

Further developments of the methods used towards a conclusive statement about 
the achievable precision of neutrino flux measurement are needed. A two-dimensional 
signal extraction procedure in the (0g,Pe) plane will allow the extraction of lepton energy 
distribution of signal events. Subsequently, a distribution of neutrino energy should 
be obtained. Development of advanced pattern recognition and fitting reconstruction 
algorithms is expected to improve on background rejection and angular and momentum 
resolutions. Finally, a dedicated analysis of systematic uncertainties should be made to 
asses the precision of the neutrino flux determination. 
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Appendix A 
Summary of cuts 



Tables A.l, A. 2 and A. 3 summarize the selection cuts for the IMD, ES~ and ES + signal 
event samples. Only the efficiencies of the selection cuts (not including reconstruction 
efficiency) are given. 



Table A.l: Selection cuts for the IMD sample. 



Cut 




Eff. after cut 


Purity after cut 


All events in detector 






0.85 x 10- 4 


Reconstruction selections 






0.19 % 


Side slabs deposit 


= MeV 


99 % 


0.25 % 


Backward deposit 


= MeV 


98 % 


0.28 % 


Vertex deposit 


< 4 MeV or 








< 500 p.e. 


98 % 


0.47 % 


(AE { ) 


< 3 MeV 


96 % 


1.4% 


max(AEi) 


< 12 MeV 


94% 


1.5% 




< 2.5 cm 


92 % 


13% 


(q/p)rec 


< 


91 % 


28 % 


Prec 


> 10 GeV 


89 % 


34% 
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Summary of cuts 
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Table A. 2: 


Selection cuts for the ES 


sample. 




Cut 


Eff. afte 


?r cut Purity after cut 


All events in detector 






2.0 x 10~ 4 


Reconstruction selections 






0.37 % 


Side slabs deposit 


= MeV 


98 % 


0.48 % 


Backward deposit 


= MeV 


95 % 


0.52 % 


Vertex deposit 


< 4 MeV or 








< 500 p.e. 


94% 


0.89 % 




< 15 cm 


89 % 


1.5% 


(q/p)rec 


< 


86 % 


7.2 % 


Prec 


< fO GeV 


73 % 


23 % 



Table A. 3: Selection cuts for the ES + sample. 



Cut 




Eff. after cut 


Purity after cut 


All events in detector 






4.1 x 10~ 4 


Reconstruction selections 






0.67 % 


Side slabs deposit 


= MeV 


98 % 


0.83 % 


Backward deposit 


= MeV 


95 % 


0.89 % 


Vertex deposit 


< 4 MeV or 








< 500 p.e. 


94% 


1.4% 


•Edev 


< 15 cm 


90 % 


2.2 % 


(<1/P)rec 


< 


87% 


24% 
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